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Abstract 

We study the quantum dynamics of relaxation, decoherence and entropy of a qubit embedded in 
an anisotropic photonic crystal (PhC) through fractional calculus. These quantum measurements 
are investigated by analytically solving the fractional Langevin equation. The qubit with frequency 
lying inside the photonic band gap (PBG) exhibits the preserving behavior of energy, coherence and 
information amount through the steady values of excited-state probability, polarization oscillation 
and von Neumann entropy. This preservation does not exist in the Markovian system with qubit 
frequency lying outside the PBG region. These accurate results are based on the appropriate 
mathematical method of fractional calculus and reasonable inference of physical phenomena. 
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I. INTRODUCTION 



A central tenet in quantum computing is to create and manipulate qubits in a controllable 
way. Many novel methods have been proposed to generate controllable entangled states 
through the surrounding environment of the qubits (lH3|- When a qubit made of a two- level 
atom is embedded inside a reservoir with structured photon density of state (DOS), the 
evolution dynamics of its states is affected by the environment and exhibits non-Markovian 
behavior flj. One o f stru c tured reservoirs , three-dhnensional (3D) photonic crystals 
(PhCs). A qubit embedded in a PhC with photon DOS containing a cut-off below the 
band edge frequency can form a photon-atom bound state and show strong non-Markovian 
dynamics which provides the basis of manipulating a qubit in this structured reservoir. 

PhCs are ordered crystals with periodic variation of refractive index and show photonic 
band gaps (PBGs) in their band structures. For a practical 3D PhC, the dispersion relation 
near the PBG edge exhibits a directional-dependent property and thus has been expressed as 
a vector form by the effective-mass approximation 7| • This 3D material with the anisotropic 
band structure possesses a DOS proportional to y/ui — cu c Q(lo — co c ), where co c is the band 
edge frequency and the Heaviside step function Q(u — u c ) characterizes the cut-off behav- 
ior. Recently, the strong effect of this structured DOS on the spontaneous emission (SE) 
dynamics is observed in the experiment of quantum dots embedded in the diamond-based 
PhC lattice |4|. The strong light-matter interaction resulting from the structured DOS leads 
to the non-Markovian SE dynamics including non-exponential decay and formation of the 
photon-atom bound state. The decay rate of this non-Markovian SE can be varied by a 
factor larger than 10 through controlling the lattice parameters of this 3D structure. This 
modification of SE rate will lead to the change of quantum coherence of the system. The 
related coherence effects include the electromagnetically induced transparency (EIT) jg], 
slow propagation of light (9], and nonlinear effects 10]. These coherent effects protect the 
embedded qubits from suffering the destruction such as entanglement sudden death (ESD) 
by the environment and make the PhCs become the promising candidate for the reservoir 
of the quantum information systems III, Il2 | 

However, the previous theoretical studies about the SE dynamics of the PhCs through the 
conventional Laplace transform method delivered the physical results with some fallacious 
predictions which are inconsistent with experiments p, Il3j . These qualitatively accurate 
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results have been used frequently as a model for studying the dynamics of quantum mea- 
surements for qubits embedded in PhCs (3), Jjj]. The misleading quantitative conclusions 
made by these studies would lead to the inaccurate strategies to manipulate qubits in PhCs. 

As a qubit is embedded in a PhC reservoir, the dynamics of this qubit is related to 
the delay Green function which is a measure of the PhC reservoir memory on the qubit. 
This delay Green function, also named memory kernel, depends strongly on the DOS of the 
PhC reservoir. The memory kernel of the PhC reservoir with a cut-off DOS exhibits the 
long-time memory effect, which leads to the non-Markovian memory dynamics of SE. This 
non-Markovian dynamics can be solved accurately through fractional calculus which have 
been applied in studying the SE dynamics of PhCs jfj, Q. The study for the isotropic PhC 
system through fractional calculus shows that the SE dynamics exhibits no photon-atom 
bound state as the qubit frequency lies outside the PBG region. This result is consistent 
with the experimental result that no prolonged lifetime effect exists for the emission peak 
lying outside the PBG region [if]]. The study about SE dynamics from anisotropic PhCs 
by the use of fractional calculus shows that the anisotropic property of the PhC reservoir 
enhances the decaying behavior of the qubit which agrees with the experimental observation 
that emission spectrum appears only in the PhC system with strong anisotropic property 
for the emission frequency lying inside the PBG region jf], \v^. These accurate results of 
the SE dynamics by the use of fractional calculus reveal that this mathematical method is 
appropriate for the PhC systems with non-Markovian dynamics. We therefore use fractional 
calculus to investigate the dynamics of quantum measurements for a qubit embedded in an 
anisotropic PhC. 

The controllable characteristic of SE dynamics from PhCs provides the basis of many 
novel methods of manipulating qubits (si, Q, 14, [lsj. For example, in Ref. [3 Bellomo et 
al. suggested high values of entanglement trapping between two independent qubits can be 
achieved through choosing the isotropic PhC materials as reservoir. This achievement of 
entanglement preservation is based on the direct link between the SE dynamics of the single 
qubit and the entanglement dynamics. Whenever the qubit forms a bound state with the 
emitted photon through the PhC reservoir, entanglement trapping follows. However, because 
they followed the predictions of previous studies with some quantitative fallaciousness joll"!^)]. 
they predicted the unphysical entanglement trapping occurred for the qubit frequency lying 
outside the PBG region. This problem also happened in Ref. [l8| where Entezar explored 
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the entanglement dynamics of a single qubit and its radiation field near the PBG edge. He 
predicted that there exists a steady-state entanglement between the qubit and radiation 
field even when the qubit frequency lies outside the PBG region. This prediction was made 
by supposing that the system would form the photon bound state as the qubit frequency 
lies outside the PBG region. This supposed fact conflicts the experimental results that 
photon bound states exist only for the region deep inside the PBG but not for the outside 



PBG region 



171 ] . Studying accurately about how a single qubit evolves in the PhC reservoir 



provides the reliable basis of strategies to manipulate qubits in PhCs. We thus aim at finding 
out the accurate dynamics of a single qubit in a PhC reservoir to provide the reliable basis 
of manipulating qubits in the PhC environment. 

In this paper, dynamics of quantum measurements including the relaxation, decoherence 
and von Neumann entropy of a qubit embedded in a PhC with one-band anisotropic model 
are investigated by the use of fractional calculus. The system exhibits distinct dynamics as 
the qubit frequency lies inside and outside PBG region. Dynamics of the non-Markovian 
system with qubit frequency lying outside PBG region exhibits quickly suppression of relax- 
ation and decoherence which implies the preservation of qubit energy and coherence. The 
coherent correlation between the qubit and the PhC reservoir leads to the preservation of 
information amount stored in the qubit through the steady entropy with non-zero value. 

The paper is organized as follows. In Sec. II, the theoretical model for the system of a 
qubit in an anisotropic PhC is depicted through the Hamiltonian and the wave function. In 
use of the fractional calculus, we express the kinetic equation of the qubit as a fractional 
Langevin equation and solve this equation analytically to construct the reduced density 
matrix of the qubit. Dynamics of relaxation, polarization, decoherence and von Neumann 
entropy are discussed based on this matrix. Quantum operations of amplitude damping, 
polarization damping and multistage computations on the single qubit in the non-Markovian 
and Markovian systems are analyzed according to these discussions. Finally, we summarize 
our results in Sec. III. 
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II. DYNAMICS OF A QUBIT IN AN ANISOTROPIC PHC 



When we consider a system consisting of a qubit embedded inside a PhC reservoir shown 
in Fig. 1, the Hamiltonian can be expressed as 

H = ifia;io(|l> (1| - |0> (0|) + + W ® 4 + /i ' c -)- (!) 

Here Co>i is the qubit transition frequency between the excited state |1) and ground state 
|0). The photon mode with frequency is created and annihilated by the operators at 
and ag. Coupling strength associated with the qubit and the photon mode with frequency 
ujj: is specified by = UJ10 ^ 10 [ 2eQ ^ v } 2 eg ■ Ud with the fixed qubit dipole moment dio = dioUd, 
sample volume V, dielectric constant eo and polarization unit vector eg of the photon mode 
with frequency wg. 

If we use the coordinate (9, 0) on the Bloch sphere to parameterize the state of the qubit 
with (0, 0) for the excited state and (n, <fi) for the ground state, then the initial state of the 
total system can be written as 

®l°£> + E^(0)|0)®|lg> (2) 

k 

with the initial coordinate (8q, 4>o) on the Bloch sphere. Here the vacuum state and one- 
excitation state of the PhC reservoir are expressed as |0g) = |0i, O2, 0g, 0g +1 , and 
|lg) = |0i, O2, •••0g_ 1 , lg, 0g +1 , ...), respectively. We assume Cg(0) = for no initial corre- 
lation between the qubit and reservoir. As the system evolves, the quantum state of the 
system in the single photon sector can be written as 

®h) + J2 c dm^\h) (3) 

k 

with initial condition u p (0) = 1, u d (0) = 1 and Cg(0) = 0. Here u p (t) and Ud(t) 
stand for the excited-state and ground-state probability amplitudes of the qubit with the 
vacuum state of the PhC reservoir while Cg(t) for the qubit in its ground state with 
one photon in the reservoir. This quantum state can also be expressed as \ip{t)) = 
j e ^(*) C0S [M]|i) + s ^[M]]o)|(g) |Og) + X)feQ(0|0) ® |lfc> with the time-dependent Bloch 
parameters [#(£), 4>{t)] when we discuss the change of the qubit's state on the Bloch sphere. 



Wo)> = 



0o 



e^cos(^)|l) + sm(^)|0) 



\m) = 



Up (t)e^cos( 9 f)\l) +u d (t)stn(^)\0) 
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By substituting the quantum state in Eq. (3) into the time-dependent Schrodinger equa- 
tion, we obtain the equations of motion for the amplitudes as 

k 

jC % {t) = 9i; u p {t)cos{e Q /2)^\ (5) 

and 

= (6) 

with detuning frequency fig = uj^ — oj\q. The last equation of motion yields Udit) = Ud(0) = 1 
meaning that the ground-state probability with vacuum photon mode will not evolve with 
time, i.e., a time constant. The other two equations can be combined as 

d f l 

-u p (t) = -J G(t- r)u p (r)dT (7) 

with the memory kernel G(t — r) = g^e~ in k^ T \ This equation of motion reveals that 
the future state of the qubit is related to the memory of the reservoir in its previous state 
through the memory kernel. As the qubit is put in free space, the memory kernel has the 
form of a Dirac delta function G(t — r) oc 8(t — r) corresponding to continuous photon 
DOS. In this case, the reservoir manifests its memory effect only at an instant time r = t 
which leads to the excited amplitude of the qubit decaying exponentially with time. This 
Markovian result in free space manifests that the qubit loses all memory of its past and 
decays quickly to its ground state. 

For the anisotropic PhC reservoir discussed here, the memory kernel manifests its memory 
effect within the entire interval (0,t) through G(t — r) = ^JJ^zii e~ t t 37r / 4 ~' 5 ( t ~ T )] with the 
coupling constant /3 1 / 2 = (uf dl Qy /uQ / (I6ire hc 3 ) and the detuning frequency 5 = cu w — co c 
of the qubit transition frequency uji from the band edge frequency u c [6(. This memory 
kernel is derived from the anisotropic dispersion relation of a practical three-dimensional 
PhC. Near the band edge frequency u c , the dispersion relation has a vector form and could 



20] ooj: m co c + A (jc — k^j where 



be expressed by the effective-mass approximation as 
the curvature A = fu c /k^ = fc 2 /u c signifies its different values in different directions 
through the scaling factor /. This anisotropic dispersion relation leads to the memory kernel 
expressed by the cut-off photon DOS p(u) through G{t-r) = J °° dco^e^^ 1 "^-^ 
with p(oj) = J-2J ^£ -^^r L ©(cJ — u c ) and the Heaviside step function <d(u — ui c ) charactering the 



cut-off behavior. With this memory kernel for the anisotropic PhC reservoir, the equation 
of motion in Eq. (7) becomes 

This equation can be further simplified through making the transformation u p (t) = e zSt U p (t) 
which gives 

d TT ^ , , rrr ^ pVW 4 [* U p (t) 



Jt u ^ +i5U ^ = ^lw^ dT - (9) 

Conventionally, this integro-differential equation is solved through Laplace transform which 
leads to the fractal phenomenon of the system from this memory kernel in Laplace image 2l| . 
This fractal phenomenon would result in the stochastic nature of the dynamical behavior of 
the system which appears as the non-Markovian dynamics. This non-Markovian dynamics 
can be solved accurately through fractional calculus which we have shown its appropriateness 
w comparing the obtained results with the experimental ones and finding their consistence 

In the following, we use fractional calculus to solve the non-Markovian dynamics of the 
system analytically. Especially, the fractional time derivative, one of the operators of frac- 
tional calculus, is used to express the integral term of the kinetic equation (9). In the 
well-known Riemann-Liouvile definition, the fractional time derivative operator d u /dt v is 



expressed as 



22|] 



d v ,/ n 1 d n f* fir) , , . 

for n — 1 < v < n and being the Gamma function. Through expressing the right- 
hand-side term of Eq. (9) as the fractional time derivative operator with order n = 1/2 
and applying the appropriate fractional operations to this kinetic equation, we arrive at the 
fractional kinetic equation of this qubit-reservoir interacting system as 

The fractional time derivative in this fractional differential equation indicates a subordinated 
stochastic process directing to a stable probability distribution 23|]. The evolution equation 
governing the future of the qubit is expressed as the form of a fractional Langevin equation 
because of its interacting with the PhC reservoir. Our aim of studying the quantum dynamics 
of relaxation, decoherence and entropy for the qubit system can be achieved by solving this 
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equation through the Laplace transform for the fractional operators. The basic formula used 
here is 

n— 1 

dt v v ' I V dt" 



Ju \ poo JU 

L \a**®\= L e- st ^f(t)dt = s»L{f(t)}-Y,s r 

m=0 



Jv—m—1 

-fit) 



dt u ~ n 



(12) 

t=o 



with s denoting the Laplace variable. The procedure of performing the Laplace transform 
on the fractional Langevin equation leads to 

U P (S) = s + iS + 2 pi/2 e i*/4 s l/2/f3/2 (13) 

where U p (s) is the Laplace transform of U p {t). This equation can be further expressed 
as a sum of partial fractions as the roots of the indicial equation Y 2 + 2/3 1//2 e" r/ ' 4 iy/ 3 / 2 + 
i8 = are found with the variable s 1 ^ 2 having been converted into Y. As /3/f 3 ^ 5, 
we have the different roots of the indicial equation which leads to the partial-fractional 
form of U p (s) = - ^ with Y x = e W< + yff^S) and Y 2 = 

e i7r/4 _ ^ _ 5 j p or ^/J 3 = 5 ) the indicial equation has degenerate root leading 

tof7 P (,) = l/(v^+^^) 2 . 

The inverse Laplace transform of these partial-fractional forms can not be found in the 
conventional mathematical table because the power of the variable s is not an integer. 
For the fractional power of s, we can find its inverse Laplace transform in the book of 
fractional calculus [24] as L" 1 j ( s J_ a ) j = Ylj=i a ^~ 1 ^t(j 1 ' ~ ^, a<i ) and L" 1 j ^\y 2 j = 

£l=i ELi a ' +m " 2 (0' + m > - 2 , a<? ) - [0" + m > - 2] ^ (0" + m > - 1> a<? )} with " = 

1,1/2,1/3,... and g = 1, 2, 1/p. Here we name the two-parameter function E t (a,a) 
as the fractional exponential function with variable t, order a and constant a. It is de- 
fined as the fractional derivative of an ordinary exponential function E t (a,a) = ^^e at = 
t a S^°=o r(i+i+i) w ^h the derivative formula ^E t (a, a) = E t (a—fM, a). The functional form 
of this fractional exponential function E t (a,a) will not be changed after being performed 
derivative operator with fractional or integral order. A linear combination of the fractional 
exponential functions is a potential solution of the fractional differential equation. By ap- 
plying these inverse Laplace transforms to the partial- fractional forms of U p (s), we obtain 
the analytical solution of the fractional differential equation [11] as 

U P (t) = 2e , ff/4 ^L_ x [^(1/2, Y 2 ) - Y?E t (l/2, Y 2 ) + Y^ - Y 2 e Y ^] , (14) 



8 



for (3 jf j£ 5; and 

03/2 i «l/2 p m/4 1 Ojffi , , tfVV*/ 4 

W«) = -^—f^^ifSI! >)-Lj—E,(\,ill// 3 ) + (l + ^)e W/'-S^-^ 

(15) 

for /3// 3 = <5. Here we have applied the recursion relation of the fractional exponential 
function E t (is, a) = aE t (u+ 1, a) + f?^qrn to these analytical expressions. The time evolution 
of the wave function is thus obtained as 

®h)+Y, c -k( t W)®\ i k) (i6) 

with Mp(t) = e lSt U p (t) expressed through Eq. (14) and (15), ua(t) = Ud(0) = 1 and 

Ed^(0| 2 = 1 " [u P (t)cos( d i)} 2 - H(y)] 2 - 

By definition, the reduced density matrix of the qubit can be directly obtained from this 
wave function through tracing over the reservoir degrees of freedom. It gives 



u p {t)e 



'*w|)|i> 



«d(0«n(y)|0> 



Pit) 



\u p (t)\ 2 cos 2 (^) i M ;(t)e-^^n(^o) 
!w p (t)e^'sm(# ) 1 - \u p (t)\ 2 cos 2 (f) 



An(*) Poo(*) 



(17) 



/ cos 2 (|i) |e^°sm(0 o )\ 
with the initial one p(0) = . The elements in this matrix 

V \e~ i(t>0 sin(0 o ) sm 2 (f) J 
are associated with the polarization and probabilities of the qubit. In the following, we 

will study the quantum measurements of relaxation, decoherence and von Neumann entropy 

based on these elements. 



A. Relaxation rate 

As the excited qubit spontaneously emits a photon, the qubit undergoes a quantum 
operation of amplitude damping, the excited-state amplitude of the qubit being reduced. 
This amplitude damping operation associated with the energy dissipation of the qubit leads 
to the flow of energy from the qubit into the environment through the emitted photon. 
How the energy flows from the initially excited qubit can be observed from the excited-state 
probability of the qubit P(t) through P(t) = \u p (t)\ 2 cos 2 (f ) = \U p {t)\ 2 (0 O = 0) shown in 



9 



Fig. 2(a). For the qubit frequency lying inside PBG [5/(3 = (coiq—co c ) / j3 < 0], the probability 
dynamics exhibits decay and oscillatory behavior before reaching a non-zero steady-state 
value. The energy flow from the qubit to the PhC reservoir is paused through this non- 
zero steady-state value. The excited-state energy of the qubit is preserved by keeping the 
emitted photon bound with the qubit. Bound states with larger steady-state values have 
less energy flow from the excited qubit to the reservoir. The rate of this energy flow can 



be measured through the relaxation rate T re i ax ,(t) = — ^fl = — jfjpr + jprrpr with U*(t) 



pn(t) 



u P (t) , f^W 



U p (t) ^ U*(t) 



being the complex conjugate of U p (t) in Eqs. (14) and (15). The dynamics of the relaxation 
rate in Fig. 2(b) from this expression decays quickly and becomes zero thereafter for the 
qubit frequency lying inside PBG (8//3 < 0). Combining with the dynamical behavior 
of excited-state probability, we find that the qubit loses partial of its energy in the very 
beginning period of time and then preserves the remaining energy in the end. The quickly 
stop of losing energy results from the formation of the bound state between the qubit and 
the reflected photon. The remaining energy is preserved through this bound state in the 
end. Energy localization around the qubit by the bound state leads to the non-decaying 
dynamics of excited-state probability in Fig. 2(a) and quickly decaying of relaxation rate 
in Fig. 2(b). Especially, as the qubit frequency is detuned deep inside the PBG region 
(6//3 = —5 and —10), the relaxation rate exhibits negative values in the initial period of 
time. This negative-relaxation-rate phenomenon reveals that the qubit can regain its energy 
in the PhC reservoir. 

On the other hand, for the qubit frequency lying outside PBG region [5/(3 = (tuio — 
oj c )//3 > 0], the relaxation rate decays slowly and remains a constant of non-zero value at 
longer time. It corresponds to that the excited qubit loses almost its energy in the initial 
period of time with (3t < 2 shown in Fig. 2(a). In this Markovian case, the emitted photon 
flees from the qubit without being re- absorbed so that it carries away all of the qubit 's 
energy to the reservoir. 

When we use a Bloch vector in the Bloch sphere to express the qubit state, the amplitude 
damping operation will perform a transformation on the components of the Bloch vector with 



the values related to the probability of losing the emitted photon (25] . For the Markovian 
system with qubit frequency lying outside the PBG region, the large-value probability of 
losing the emitted photon will convert the Bloch vector toward the ground state point on the 
Bloch sphere. The small- value probability of losing the emitted photon in the non-Markovian 



10 



system with qubit frequency lying inside PBG leads to the Bloch vector towards the non-zero 
large surface area near the ground-state point on the Bloch sphere. The contracting effect 
of the amplitude damping on the Bloch vector leading to the contracted surface area of the 
Bloch sphere available for the Bloch vector is suppressed in the non-Markovian system. 



B. Decoherence 



As a qubit interacts with the environment, there will randomize the polarization of 
each qubit. The quantum phase information of a qubit carried by the qubit polariza- 
tion will thus escape from the qubit into the environment through this randomization of 
polarization and lead to the quantum decoherence. As the qubit undergoes the polariza- 
tion damping operation through this quantum decoherence, the off-diagonal elements of 
the reduced density matrix will decay with time. The completely losing of the quantum 
coherence of one qubit will lead to the occurrence of ESD in the multiple-qubit system. 
That is, the time period in which the entanglement could be usefully exploited in the 
multi-qubit system is limited. This limitation is related to the dynamics of qubit polar- 
ization and decoherence. Here we show these dynamics in Fig. 3 through the expres- 
sions of qubit polarization Pz(t) = \\pwit) + Poi(t)} = Re[e l5t U p (t)] and decoherence rate 



u P (t)+u;(t) 



pio(t)+p i(t) [u p (t)+u;(t) 
In Fig. 3(a), the polarization dynamics of the Markovian system with qubit frequency 

lying outside the PBG region (S//3 = 2) exhibits fast damping behavior which implies that 
the qubit loses all its polarization in the initial period of time with (3t < 4. This polarization 
damping results from the elastic scattering by the incoherence photons of the environment. 
The elastic scattering randomizes the orientation of the qubit's polarization and leads to 
the decay of the off-diagonal elements of the density matrix. The fast damping of the qubit 
polarization reveals that the probability of scattering the qubit by the incoherent photons 
in the reservoir is large in the Markovian system. The continuous photon DOS outside the 
PBG region [see Fig. 1(c)] provides the free-space-like environment to scatter the qubit. 

In contrast, the polarization dynamics of the non-Markovian system with qubit frequency 
lying inside PBG (5/ j3 < 0) exhibits non-decaying oscillation. The qubit loses partial of its 
polarization in the very beginning period of time and then preserves the remaining polar- 
ization through the steady oscillation. This steady polarization ensures that the probability 
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of scattering the qubit by the incoherent photons is greatly lowered in the non-Markovian 
system. The inhibited photon DOS of the PhC reservoir inside the PBG region [see Fig. 
1(c)] is expected to lower this probability of scattering the qubit and lead to the preservation 
of phase information of the qubit through the steady polarization. States with larger ampli- 
tudes of steady oscillation have less lose of phase information from the qubit to the reservoir. 
The rate of this lose of phase information can be measured by the decoherence rate Vdec^t) 
shown in Fig. 3(b). The decoherence rate of the non-Markovian system decays quickly to 
zero. Combining the results of Fig. 3(a), we infer that the phase information in the qubit is 
partially lost in the very beginning period of time and then is preserved through the steady 
polarization. In this non-Markovian case, the interaction with the PhC reservoir does not 
result in the total lose of the coherence but in the preservation of the partial coherence 
instead. This coherence preservation reveals that PhC materials can provide a promising 
environment for the purpose of information processing because they meet the requirement 
of engineering the information-processing systems through the greatly minimized quantum 
decoherence. 

Comparing the dynamics of the relaxation rate [Fig. 2(b)] with that of the decoherence 
rate [Fig. 3(b)], we can find the similarity between these two quantum measurements. 
This is because both the relaxation and decoherence arise from the effects of spontaneous 
emission. As the emitted photon flees from the qubit, it carries away both the energy and 
phase information of the qubit. We compare these two quantum measurements through 
calculating the ratio of the decoherence rate to the relaxation rate. This ratio exhibits 
constant at 0.5 for both the Markovian and non-Markovian regimes. This result shows that 
the decoherence will proceed much slower than the energy relaxation in this PhC reservoir. 
That is, the decoherence time 7\ is twice the time of the energy relaxation T2, i.e., T\ = 2T2. 
This familiar relation between the decoherence time and relaxation time for the generic case 
of a two-level system was also observed in a transverse reservoir with a Lorentzian spectrum 
in the infinite-temperature limit 25J, |26[ . When the energy relaxation is completed, almost 
half of the initial coherence of the qubit will still survive. 
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C. Von Neumann entropy 



Entropy, a key concept of quantum information theory, measures how much informa- 
tion exists in the state of a physical system. The amount of information will be changed 
as the state is correlated to the environment. The correlation between the environment 
and the state will transform the initially pure state into a finally mixed state where 
the amount of information of the state is changed. We thus can estimate the correla- 
tion between the qubit and reservoir through the entropy of the quantum system, i.e., 
von Neumann entropy. Von Neumann entropy for a quantum state p(t) is defined as 
S(t) = —Tr[p(t)logp(t)\ = -^^'05^ with A, being the eigenvalues of the matrix 
p(t). In this expression, the entropy quantifies the optimal compression for the physical 
resources required to solve some information problem because in the measurement the 
eigenvalues \ stand for the expected values and log\ for the probability of obtaining 
Aj. That is, von Neumann entropy predicts the upper bound of the information we will 
gain after we measure a quantum state. This gain of information after measurement cor- 
responds to the information amount of the state before measurement. The value of this 
entropy stands for the degree of correlation between the qubit and the PhC reservoir. It 
reaches the maximal value S max = logd, d being the dimension of the matrix p(t), for the 
qubit's maximally mixed state and zero value for the pure state. With the eigenvalues 



\± = \ ll± Jl- 4cos 4 (f ) [\u p (t)\ 2 - \u p {t)\ 4 ] \ of the matrix in Eq. (17), we show the 



von Neumann entropy in Fig. 4 for the initially excited qubit (6> = 0). The entropy has 
its minimal value zero at t — and reaches its maximal value log2 = 0.693 at the very 
beginning of time. After a period of time on the order of the decay timescale, the entropy 
becomes steady with nonzero value in the non-Markovian system (5/(3 < 0) and zero in the 
Markovian system (5/(3 = 2). These results show that the initially pure system becomes 
maximally mixed in the very beginning period of time. As the qubit equilibrates with the 
PhC reservoir, the system becomes less mixed. In the Markovian system, the qubit recovers 
its initially pure state as it quickly equilibrates with the reservoir. However, in the non- 
Markovian system, the coherent correlation with the PhC reservoir leads to the preservation 
of the mixed state. The greatly lowered probability of the qubit scattering by the incoherent 
photons in the PBG reservoir preserves the coherence and information amount of the qubit. 
When we perform a multi-stage computations on the available information through the 
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Markovian chain, the data processing inequality of von Neumann entropy, a basic inequality 
of information theory, states that the information about the output of the computation will 
decrease with time if the qubit storing the information is correlated to the Markovian en- 



vironment 



25j . This statement agrees with the result we obtained here. In the Markovian 



system, the newly produced mixed state by the correlation between the qubit and reservoir 
recovers its initially pure state and loses the information stored in the qubit. The informa- 
tions yielded by the multi-stage computations will thus be independent of each other. On 
the other hand, in the non-Markovian system, the newly produced mixed state will preserve 
the information about its output of the previous stage. The memory effect of the PhC 
reservoir on the qubit 's previous state leads to the preservation of the stored information in 
the mixed-state qubit. Further computation operations on this mixed state can be used to 
increase the amount of mutual information between the outputs of the operations and the 
previous-stage information about the qubit state. 



III. CONCLUSION 



We have studied the dynamics of quantum measurements of relaxation, decoherence and 
von Neumann entropy for a qubit in an anisotropic PhC. Through applying fractional cal- 
culus to solve the fractional Langevin equation, we obtain the analytical expressions of 
these quantum measurements. The quantum operations corresponding to these quantum 
measurements act differently in the non-Markovian and Markovian systems which are dif- 
ferentiated according to the qubit frequency lying inside or outside the PBG region. For the 
non-Markovian system with qubit frequency lying inside PBG (S//3 < 0), the relaxation dy- 
namics shows that the qubit can preserve the remaining energy through forming the bound 
state with the photon. The small-value probability of losing the emitted photon leads to 
the suppressed contracting effect of the amplitude damping operator on the qubit 's vector 
in the Bloch sphere. Without forming the bound state with the emitted photon, the excited 
qubit of the Markovian system with frequency lying outside the PBG region (6 /ft > 0) loses 
all its energy to the reservoir. This large-value probability of losing the emitted photon 
leads to the Bloch vector of the qubit contracting to the ground-state point on the Bloch 
sphere. The distinct behavior between the non-Markovian and Markovian systems also ex- 
ists in the decoherence and entropy dynamics. Steady polarization oscillation and entropy 
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with non-zero value in the non-Markovian system reveals that the greatly lowered proba- 
bility of scattering the qubit by the incoherent photons in the PhC reservoir leads to the 
preservation of the coherence and information amount of the qubit. This preservation does 
not exist in the Markovian system with continuous photon DOS outside the PBG region 
where the PhC reservoir provides the free-space-like environment with large probability of 
scattering the qubit. Scattering by the incoherent photons leads to fast damping of qubit 
polarization and losing the information stored in the qubit. These different results from 
those of previous studies are based on the appropriate mathematical method of fractional 
calculus and reasonable inference of physical phenomena. 
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FIG. 1: (Color online) (a) A qubit with excited state |1) and ground state |0). The transition 
frequency wio is nearly resonant with the frequency range of the PhC reservoir, (b) Directional 
dependent dispersion relation near band edge expressed by the effective-mass approximation with 
the edge frequency uj c . (c) Photon DOS p{u) of the anisotropic PhC reservoir exhibiting cut-off 
photon mode below the edge frequency u c . 
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FIG. 2: (Color online) Dynamics of (a) the qubit's excited-state probability P(t) and (b) relaxation 
rate T re i ax .{t) of the qubit with different detuning frequencies S//3 = (wio — 0J c )/P from the band 
edge frequency co c of the PhC reservoir. 
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FIG. 3: (Color online) (a) Polarization Pz(t) = ^[pio(t) + poi(t)] and (b) Decoherence rate 
r<fe c .(i) = - ^°jij+^(tj of the q ubit with frequency lying inside (S//3 < 0) and outside (6/0 = 2) 
the PBG region . 
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FIG. 4: (Color online) Von Neumann entropy S(t) = —Tr [p(t)log p(t)] 
the non-Markovian (6/(3 < 0) and Markovian (5/(3 = 2) systems. 
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